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PRELIMINARY  RADAR  FEATURE  EXTRACTION 

AND 

RECOGNITION  USING  TEXTURE  MEASUREMENT 


INTRODUCTION 

Recently  an  image-domain  processing  technique  was  investigated  and  implemented 
with  an  experimental  solid-state  sensor  array-minicomputer  system  at  the  U.S.  Army 
Engineer  Topographic  Laboratories  (ETL).1  The  system  employs  a  32-element  by  32- 
element  solid-state  sensor  array  to  convert  images  into  electronic  signals  and  a  minicom¬ 
puter  to  process  the  signals  for  extracting  and  classifying  the  cartographic  features  from 
the  imagery  into  preassigned  categories  based  on  a  feature  vector.  The  images  under 
test  and  investigation  were  selected  aerial  photographs. 

The  purpose  of  this  effort  was  to  modify  and  verify  the  above  system  for  ex¬ 
tracting  and  recognizing  features  from  a  selected  set  of  radar  imagery  of  the  Huntsville. 
Alabama,  area.  Description  of  the  system  modification  is  followed  by  a  discussion  of  the 
selection  of  feature  vector  components  and  classification  strategy.  Classification  results 
for  a  set  of  selected  radar  imagery  are  presented.  Finally,  conclusions  are  given  together 
with  comments  regarding  extensions  of  this  work. 


SYSTEM  DESCRIPTION 

The  hardware  portion  of  the  system  used  for  this  experimentation  is  essentially 
the  same  as  the  one  reported  previously.2  The  voltage  of  the  light  source  was  increased 
because  most  of  the  radar  imagery  was  rather  dark.  A  new  software  program  was  de¬ 
veloped  for  the  extraction  of  radar  features  because  radar  signatures  of  terrain  features 
are  quite  different  from  their  counterparts,  the  cartographic  features  from  aerial  photo¬ 
graphs. 


The  block  diagram  of  the  system  is  shown  in  figure  1 .  A  9-inch  by  9-inch  glass 
plate  mounted  with  strips  of  radar  imagery  is  illuminated  by  a  white  light  source,  and 
a  section  of  the  image  is  projected  onto  a  Reticon  32-element  by  32-element,  solid- 
state  array  through  an  imaging  lens.  The  array  converts  the  optical  energy  of  the  image 
into  a  video  signal.  The  video  signal  is  quantized  into  10  bits  of  digital  signals  and  sent 
to  the  Hewlett-Packard  2108  minicomputer  for  processing.  The  computer  first  takes  in 
the  quantized  signals  of  32  pixels  by  32  pixels  of  1,024  gray  levels  array.  With  the 
brightest  and  darkest  pixels  within  a  frame  as  the  maximum  and  minimum,  this  quantized 
image  array  is  next  scaled  down  to  become  32  pixels  by  32  pixels  of  16  gray  levels. 
The  joint  probability  matrix  of  this  scale  array  is  then  obtained.  The  next  step  is  to  com¬ 
pute  a  feature  vector  based  on  the  joint  probability  matrix.  Although  nine  feature  vector 
components  were  computed,  only  two  are  needed  for  classification  of  the  selected  radar 
imagery.  A  sequential  template-matching  classifier  is  used  to  classify  the  input  images 
into  one  of  the  preassigned  image  categories;  an  image  is  recognized  as  a  reject  if  it  does 
not  belong  to  any  of  these  categories.  The  classification  result  is  then  indicated  on  a  CRT 
console.  At  the  end  of  classification  a  signal  is  sent  to  the  translational  stage  controllers 
to  move  the  stages  in  the  predetermined  X  and  Y  positions,  and  a  new  section  of  image 
is  projected  onto  the  surface  of  the  solid-state  array.  The  procedure  described  repeats 
until  all  preselected  image  sections  are  classified. 


RADAR  FEATURE  EXTRACTION 


The  system  described  in  the  previous  report  uses  image  histogram  and  image 
texture  as  the  bases  for  developing  a  feature  vector.3  Only  the  image-texture  technique 
was  considered  for  the  extraction  of  radar  features  from  a  selected  set  of  radar  imagery 
because  the  image  categories  of  interest  are  easier  to  separate.4 


Image-texture  features  (second-order  image  statistics)  are  based  on  the  definition 
of  the  joint  probability  distribution  of  pairs  of  pixels.  Pratt  stated  that  the  two-dimen¬ 
sional  histogram  can  be  considered  as  an  estimate  of  joint  probability  distribution.5 
Consider  a  pair  of  pixels  F(j,  k)  and  F(m,  n)  that  are  separated  by  7  radial  units,  and 
are  at  an  angle  6  with  respect  to  the  x-axis  of  the  measurement  window.  The  histogram 
estimate  of  the  second-order  distribution  is  given  by  Pratt6  as 


P(a.  b) 


N(a.  b) 
M 


(1) 


Where  M  is  the  total  number  of  all  occurrences  in  the  measurement  window  and  N(a.  b) 
denotes  the  number  of  occurrences  for  which  FO'.  k)  =  a.  F(m,  n)  =  b.  Various  texture 
measures  that  have  been  used  in  this  study  are  listed  in  appendix  A  (and  also  in  ETL-- 
0297). 7  These  measures  are  as  follows: 


1. 

Mean 

2. 

Variance 

3. 

Covariance 

4. 

Autocorrelation 

5. 

Absolute  Value 

6. 

Energy 

7. 

Inverse  Difference 

8. 

Inertia 

9. 

Entropy 

3P.I  .  Chon,  A  Sensing  Array  System  with  Image  Statistics  Processing.  U.S.  Army  1-ngineer  Topographic  Laboratories, 
f  ort  Bclvoir,  Virginia.  LTL-0297.  May  1982,  AD-A119  259. 

4W.K.  Pratt.  Digital  Image  Processing.  Now  York,  John  Wiley  and  Sons,  Inc..  1978. 

SIhiJ. 

6lhid. 

7Chen,  op.  t  it. 


For  our  application.  P(a.  b)  was  made  to  be  symmetrical  so  that  a  =  T>.  and 
Va  =  Vfe  (see  appendix  A).  Each  input  image  was  first  scaled  down  from  1,024  to  16 
gray  levels  (L  =  16).  Nine  components  of  the  feature  vector  based  on  these  equations 
given  in  appendix  A  and  equation  (1)  were  then  computed.  Equation  (1)  was  eva¬ 
luated  for  0  values  of  0.45.  90,  and  135  degrees.  The  corresponding  feature  vector 
components  of  different  O' s  for  each  image  category  of  interest  were  compared.  It 
was  discovered  that  only  two  components  of  the  feature  vector,  namely  the  covariance 
and  the  autocorrelation,  and  the  number  of  pixels  within  the  measurement  window  that 
are  above  a  specific  threshold  value  (NPATV)  were  required  for  classification  purposes. 


CLASSIFIER 

For  the  selected  set  of  radar  imagery,  only  two  components  of  the  feature  vecto 
computed  in  the  previous  section  plus  the  number  of  pixels  above  a  threshold  value 
(NPATV)  were  used  to  constitute  a  three-dimensional  sequential  template-matching 
classifier.  These  three  components  are  as  follows: 

1 .  Covariance 

2.  Autocorrelation 

3.  The  number  of  pixels  above  the 
threshold  value  (NPATV). 

Many  prototype  image  samples  were  obtained  from  a  set  of  radar  imagery  to  de¬ 
termine  the  upper  and  lower  limits  of  the  template  values  for  each  image  category.  Three 
template  ranges  for  the  covariance  were  defined  as  follows:  water.  -2.0  to  1.5;  field. 
0  to  1.5;  and  city  and  forest.  1.5  to  6.  The  template  ranges  for  the  autocorrelation  were 
designated  to  be  city.  0  to  40  and  forest,  above  40  to  150.  Likewise,  the  template  ranges 
for  the  NPATV  were  set  to  be  water.  0  to  500  and  field,  above  501  to  1023  (see  figure 
2). 

The  covariance  of  the  unknown  incoming  input  image  is  first  compared  to  the 
template  ranges  of  the  covariance  template.  If  its  value  is  within  the  range  of  -2.0  to 
1.5.  then  the  NPATV  of  the  input  image  is  compared  to  its  corresponding  template. 
If  it  is  equal  to  or  less  than  500.  the  input  image  is  classified  as  “water.”  If  NPATV 
is  greater  than  500  or  the  covariance  is  not  within  the  range  of  -2.0  to  1.5.  then  the 
next  test  will  be  performed.  The  sequence  of  tests  is  always  from  “water”  to  “field” 
to  “forest”  and  finally  to  “city”  as  listed  in  appendix  B.  The  covariance  template  is 
used  as  the  preliminary  template  for  all  four  image  categories  of  interest.  The  auto¬ 
correlation  template  is  selected  as  the  final  decision  template  for  “forest”  and  “city.” 
while  the  NPATV  template  is  employed  as  the  final  decision  template  for  "water" 
and  “field.”  Finally,  if  the  unknown  image  does  not  belong  to  any  step  of  the  test 
described,  it  is  then  classified  as  “not  recognized.” 


FIGURE  2.  Template  Ranges  for  Feature  Vector  Used  for  Radar  Image  Classification. 


TEST  RESULTS 


A  set  of  high-quality,  scale  of  1  to  100,000.  X-band,  synthetic  aperture  radar 
imagery  from  the  Huntsville,  Alabama,  area  was  used  for  this  experimentation.  The  set 
consisted  of  image  categories  such  as  city  (combination  of  commercial  and  residential 
structures,  DLMS  categories  #504  F1C  301  and  #505  F1C  401).  field  (agriculture, 
used  primarily  for  crop  and  pasture  land.  DLMS  category  #510  F1C  950).  water  (river, 
smooth  fresh  water.  DLMS  category  #501  F1C941).  and  forest  (deciduous.  DLMS 
category  #510  FIC  952).  Figures  3(a)  to  6(a)  show  the  line  printer  output  of  the 
typical  radar  image  categories  of  city,  field,  water,  and  forest  respectively.  Each  is  printed 
in  16  gray  shades  by  line  printing.  In  figures  3(b)  to  6(b),  the  nine  feature  vector 
components  computed  by  using  equations  in  appendix  A  and  equation  ( I )  and  the 
NPATV  are  shown  for  each  image  category. 

Approximately  100  images  covering  all  four  categories  were  scanned.  These  images 
were  used  as  input  to  evaluate  the  classification  accuracy  of  the  sequential  template¬ 
matching  classifier.  The  result  is  illustrated  in  figure  7.  An  overall  classification  accuracy 
of  approximately  92  percent  was  obtained.  The  category  determination  of  prototype 
(or  reference)  images  that  were  used  for  this  classifier  was  based  on  the  ground  truth 
located  from  a  map  of  the  same  area. 

The  variations  of  the  feature  vector  components  with  respect  to  the  texture-mea¬ 
surement  angle,  0,  for  the  selected  radar-image  categories  of  city,  field,  water,  and 
forest  are  illustrated  in  figures  8  through  11.  respectively.  The  covariance  measure 
for  city  is  highly  directional,  as  indicated  in  figure  8.  The  covariance  for  the  horizontal 
and  vertical  (0  =  0  and  90  degrees)  directions  is  almost  twice  that  measured  in  the 
diagonal  sense  (0  =  45  and  135  degrees).  The  big  jump  in  covariance  of  water  for  the 
measurement  direction,  6  =  90  degrees,  was  due  to  a  thin,  long  bright  object  running 
through  in  the  vertical  direction  in  that  particular  measurement  window.  Other  variations 
are  relatively  small  and  insignificant.  The  classification  accuracy  shown  in  figure  7  was 
obtained  by  using  the  texture-measurement  angle.  0  =  0  degree. 
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FIGURE  3.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  Components  for  City. 
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FIGURE  4.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  Components  for  Field. 
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FIGURE  5.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  Components  for  Water. 
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FIGURE  6.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  Components  for  Forest. 
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FIGURE  7.  Classification  Results  for  Sequential  Template  Matching  Classifier. 
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INVERSE  SECOND  ORDER  MOMENT  *  .4.5073 

ENTROPY-  3.413*1 
CONTRAST-  6.52029 
ABSULUTE  VALUE-  1.83767 
AUTOCORRELATION-  31.49950 
NO.  UP  ELEMENTS) IT-  815 
SECOND  ORDER  STATISTICS 
NFP  T'JL  1  1981  IS:  26 

CTTY<  0-9*1  DECREES) 

SUM  PROBABILITY-  1.09* 

AVE-  5.3V4  VAR—  S . ' 64  COV-  3.08291 
ANGULAR  SECOND  ORDER  MOMENT  •  (14282 

INVERSE  SECONU  ORDER  MOMENT  -  .52379 

ENTROPY-  3.71*20 
PONTRAST-  4.1*129 
ABSOLUTE  VALUE-  1 . 40527 
AUTOCORRELATION-  31.9*271 
NO.  OP  ELEMENTS) IT-  814 
SECOND  ORDER  STATISTICS 
4PP  HJL  1  1981  l  S :  27 

CTTY< 0-135  DEGREES) 

SUM  PROBABILITY-  1  .  0(10 

AUE*  5.454  VAR-  4.929  COV-  1.77023 
ANGULAR  SECOND  ORDER  MOMENT  >  .  (1350* 

INVERSE  SECOND  ORDER  MOMENT  -  44541 

ENTROPY-  3.88526 
CONTRAST-  6.31842 
ABSOLUTE  VALUE-  l .  796116 
AUTOCORRELATION-  31.51303 
NO.  OP  ELEMENTS) IT-  819 
SECOND  ORDER  STATISTICS 


FIGURE  8.  Variation  of  Feature  Vector  Components  with  Respect  to  0  for  City. 
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UEO  IIJL  t  1981  IS 1 28 
PTELD<0-0  DEGREE) 

SUM  PROBABILITY-  1 . 0U0 

AVE»  4.M?  VAR-  1.443  CQV-  .  96662 
ANGULAR  SECOND  ORDER  MOMENT  -  .3123S 

INVERSE  SECOND  ORDER  MOMENT  -  .78491 

ENTROPY-  t. 94701 
CONTRAST a  . 98262 
ABSOLUTE  VALUE ■  .50706 

AUTOCORRELATION-  22.7*897 
NO.  OP  ELEMENTS > IT-  139 
SECOND  ORDER  STATISTICS 
OPT)  TIJL  1  I98t  15:29 

PTFLD<  i)«4S  DECREES  ) 

SUM  PROBABILITY-  i.1'09 

AV*«  4.608  VAR ■  1 . 4S0  COU-  . 69209 

ANGULAR  SECOND  ORDER  MOMENT  »  30S67 

INVERSE  SECOND  OROER  MOMENT  -  .76567 

ENTROPYa  3. 01 ISO 
CONTRASTa  1,51509 
ABSOLUTE  VALUE-  .60146 
AUTOCQRRELA riON-  21 . 92299 
NO.  OP  ELEMENTS) IT-  13’ 

SECOND  ORDER  STATISTICS 
ypn  JIJL  l  1981  15-.29 

PTELD<  0-90  DEGREES) 

SUN  PROBABILITY-  l . 000 

AVEa  4.614  VAR-  1.619  COVa  1.13963 
ANGULAR  SECONO  ORDER  MOMENT  a  .20983 

INVERSE  SECOND  ORDER  MOMENT  -  .80782 

ENTROPY-  1 . 98769 
CONTRASTa  . 98367 
ABSOLUTE  VALUE-  .45867 
AUTOCORRELATION-  22 . 4324- 
NO  .  OP  ELEMENTS) IT-  139 
SECOND  ORDER  STATISTICS 
UED  JUL  l  1981  15:30 

FTPLU< 9-135  DEGREES) 

SUM  PROBABILITY-  1.000 

AVE*  4.589  VAB-  1.447  COV-  .49994 
ANGULAR  SECOND  ORDER  MOMENT  •  .20528 

INVERSE  SECOND  ORDER  MOMENT  -  . 7S907 

ENTROPY-  1.99932 
CONTRAST-  1 . 49428 
ABSOLUTE  VALUE-  .61186 
AUTOCORRELATION-  21 . 75858 
NO.  OP  ELEMENTS) IT-  142 
SECOND  OROER  STATISTICS 


FIGURE  9.  Variation  of  Feature  Vector  Components  with  Respect  to  0  for  Field. 
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MED  HJL  1  1981  15:  ta 

MATES <4-4  DECS EE ) 
siim  b40Bability«  1.444 

AVE-  11.411  UAB-  4  94(1  COV-  .39838 
ANCULA*  SECOND  OBOES  MOMENT  -  . 1 1264 

t NOES SC  SECOND  OBOES  MOHEHT  -  .51294 

ENTROPY-  2.66173 
CONTRAST-  8.(14383 
ABSOLUTE  VALUE*  1.6*496 
AtJTOCOB  B  EL  A  T I  ON-  135  .  72278 
NO.  OF  ELEMENTS) TT-  0 

..SECOND  OBOES  STATISTICS 
'mFD  riJL  1  1981  15:  19 

MATER <4-45  DECS  EES ) 

sum  probability-  i.on# 

AWE-  12.240  VA*-  5.371  COW-  (.23459 
ANCULAS  SECOND  OBOES  MOMENT  a  .18252 

INVERSE  SECOND  OBDER  MOMENT  -  .54133 

ENTBOPY-  2.73039 
CONTRAST-  8 . 27367 
ABSOLUTE  VALUE-  1 . 7*796 
AUTOCORRELATION-  154 . 4822« 

NO.  OF  ELEMENTS) fT«  4 

SECOND  OBOES  STArlSTiee 
MED  IIJL  l  1981  1 5 1  l  - 

NATES  <  4—98  OECSEES) 

SUM  PROBABILITY-  1.0114 

AWE-  14.641  WAS-  6.124  COW-  3.51939 
ANCULAS  SECOND  OBOES  MOMENT  -  .11907 

INWEBSE  SECOND  OB DEB  MOMENT  -  55414 

ENTBOPY-  .7.74480 
C0NTBA8T-  5.20948 
ABSOLUTE  VALUE-  1 . 40776 
A(  ITOCOBBELAT ION-  114.75345 
NO.  OF  ELEMENTS) IT-  4 

5FCONO  OBDEB  STATISTICS 
MET)  IIJL  l  (981  15:24 

MATER < 4-135  DECREES > 

SUM  PROBABILITY-  1.908 

AWE-  11.523  VAP-  4.746  COW-  85443 
ANCULAS  SECOND  OP DEB  MOMENT  •  .19713 

INVERSE  SECOND  OBDEB  MOMENT  -  .47419 

ENTBOPY-  2.67197 
CONTRAST-  7.79*47 
ABSOLUTE  VALUE-  1.77944 
A4JTOCOPBELATTON-  133.64319 
NO.  UP  ELEMENTS) IT-  0 

5ECUN0  ORDER  SVATISTirs 


FIGURE  10.  Variation  of  Feature  Vector  Components  with  Respect  to  0  for  Water. 
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Wffn  TIJL  l  1991  15:  22 
FO*EST<*-*  DEGREE) 
sum  probability-  1.900 

AVC-  7.0  A*  UM*  3.59S  COW-  2.3335A 
ANGULAR  3ECUND  ORDER  MOMENT  •  .  0353 A 

INVERSE  SECOND  Off  Mill  MOMENT  -  .55526 

ENTRANT-  3.67226 
CONTRAST-  2. 5*2*2 
ABSOLUTE  VALUE-  1.13306 
AUTOCORRELATION*  51.9949* 

NO.  Off  ELEMENTS > IT-  Y66 

_5FC0ND  ORDER  STATISTICS 
wffO  NJL  l  1951  15:23 

FOREST < H— AS  DEGREES) 

3UM  PROBABILITY-  1.0M0 

AWE-  7.051  WAR-  1 . 447  COW-  1.84940 
ANGULAR  SECOND  ORDER  MOMENT  -  . 92932 

INVERSE  SECOND  ORDER  MOMENT  -  .47912 

ENTROPY-  3.9646* 

CONTRAST-  3 . 59*25 
ABSOLUTE  UALUE-  1 . 4339? 

AUTOCORRELATION-  SI. 56502 
NO.  OF  ELEMENTS) IT-  946 
.5ECQND  ORDER  STATISTICS 
OFI)  riJL  1  .9*1  15:24 

FOREST ( <l—9<!  DECREES) 

SUM  PROBABILITY-  1.90# 

AWE-  6.953  WAR-  4.224  COW-  2.7*104 
ANGULAR  SECOND  ORDER  MOMENT  -  .02*94 

INVERSE  SECOND  OkOCR  MOMENT  -  .51132 

ENTROPY-  3.91*64 
CONTRAST-  2. 96S93 
ABSOLUTE  VALUE-  1.2731* 

AUTOCORRELATION*  51 . 55B49 
NO.  Off  ELEMENTS > IT*  966 
SECOND  ORDER  STATISTICS 
OF!)  TUL  l  1*91  15 1 24 

smresti 9-135  DEGREES) 

SUN  PROBABILITY—  1.10* 

AVf-  7. *72  WAff-  3.439  COV-  1.597*9 
ANGULAR  SECOND  ORDER  MOMENT  -  .*3779 

INVERSE  SECOND  ORDER  MOMENT  -  .44329 

CNTROffY—  3.99*63 
CONTRAST-  4.19199 
ABSOLUTE  VALUE-  1.52341 
AUTOCORRELATION-  51.59937 
NO.  Off  ELEMENTS) IT-  945 
SFCONO  ORDER  STATISTICS 


FIGURE  11.  Variation  of  Feature  Vector  Components  with  Respect  to  0  for  Forest. 
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CONCLUSIONS 


1.  The  image-texture  technique  provides  an  effective  means  for  evaluating  texture 
and  coarseness  of  radar  area  features. 

2.  The  technique  is  most  applicable  for  extracting  and  classifying  if  the  search  win¬ 
dow  contains  only  a  single  category  of  radar  features  such  as  city,  forest,  water, 
or  field.  Multiple  categories  of  radar  features  contained  in  a  search  window  were 
mostly  misclassified  or  rejected  as  not  recognized.  Determination  and  detection  of 
boundaries  between  different  radar  features  are  subjects  of  future  research. 

3.  A  preliminary  classification  accuracy  of  slightly  better  than  90  percent  was 
obtained  for  a  selected  set  of  radar  imagery  from  the  Huntsville.  Alabama,  area. 

4.  The  technique  will  be  extended,  and  similar  experiments  will  be  conducted  for  a 
wide  range  of  radar  imagery  from  various  locations  and  for  imagery  taken  with 
different  radar  depression  angles. 
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APPENDIX  A.  Feature  Vector  Components 

L-l  L  - 1 

Mean:  a  =  £  ]T  aPfa.b) 

a  =  0  b  =  0 


L-l  L-l 

b  =  Y.  Z  bPfa.b) 

a  =  0  b  =  0 

L-1  L-l 

Variance:  Va  =  £  £  (a  -  a)2  P(a,  b) 

a  =  0  b  =  0 


L-l  L-l 

vb  =  Z  Z  (b-b)2  P(a,b) 

a=0  b=0 


L-i  L-l 

Covariance:  CQ  =  £  Z  (a -a)  (b-b)  P(a,  b) 

a=0  a=0 


L-l  L-l 

Autocorrelation:  Au  =  £  a  b  P(a,  b) 

a=  0  b  =  0 


L-l  L-l 

Absolute  Value:  Ab  =  £  ^  la-b  I  P(a,  b) 

a=0  b=0 


L-l  L-l 

Energy:  Eg  =  £  £  [P(a,  b)]2 

a  =  0  b  =  0 


APPENDIX  A.  Feature  Vector  Components  (Continued) 


Inverse  Difference: 

Id 

L  —  1 

=  e 

L  —  1 

E 

P(a,  b) 

a  =  0 

b  =  0 

+  (a  -  b)2 

L  —1 

L  —1 

Inertia: 

In 

=  E 

E 

(a  -  b)2  P(a,  b) 

a  =  0 

b  =  0 

L  —l 

L  — l 

> 

Entropy: 

En 

=  ~E 

E 

T(a,  b)  Iog2  f  Pf a,  b)J 

a  =  0 

a  =  0 
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APPENDIX  B.  Computer  Printout 


S:  F  V  C  T  3  T  = '}  O  0  0  4  IS  OH  C  R  0  0  •>  0  9  USING  0000b  8  L  K  S  R  =  0  0  0  0 

v v Vi  ?TH4  j  L 
0  0  0  2  C 

0003  c************SUBRQUTINE  '  F  D  C  T  3  ‘  - -  REV  07/25/80**************** 

0  0  0  4  C 

0005  C  SUBROUTINE  TO  PERFORM  FEATURE  CLASSIFICATION  FOR  PROGRAM 

•)  0  0  i  c  ••  R  D  T  X  1  - 

0  0  7  C 

0  0  0  S  SUBROUTINE  FDCT3( CO  V, AU ,  ITS . LUOT > 

0  00  9  IF(COV  LT.-2.0R.C0V.GT.1  ,50)GU  TO  810 

0  0  i  0  I  F  <:  ITS  IT. 0.  OR. ITS. GT. 500)  GO  Tu  810 

0  0  i  i  y  R  I  1 1  (  LUOT,  8 > 

v  0  12  S  0  I  0  S  8  8 

ooi3  810  IFCCQV  LT.O.OR.COV.GT.l  5  0  )  G 0  TO  820 

0014  I F<  I TS . L T  500  .  OR  .  I TS  .  GT  .  1024  >G0  TO  320 

0015  UR  I  TEC  LUOT, 81  ) 

001b  GO  TO  S88 

0  0  1  ?  820  IFCCGV.LT  1  50  .  OR  .  COV  GT  .  fc  0  >00  TO  330 

0  018  IFlAU.LT.  40.0R.AU  GT  2  0  0  >  G  0  TO  33  0 

0019  UR  I  TEc.  LUOT,  82  ) 

0020  GO  TO  883 

002  1  830  I  r  <  CuV  LT  1  .50 .OR. COV  GT .8  ) G 0  TO  840 

■  ‘  0  2  2  !F«:aU  LT  .  00  .OR  .  A  U  .  GT.  40  )G0  TO  840 

0  0  2  3  U  S I T  E ( L  U  0  T / 3  3  ) 

•> 2  4  GO  TO  88  3 

0  02  5  8>»o  WRITE'-  LUOT,  84  ) 

CM.  2  8  GO  TO  88  8 

0  02  7  90  F  0  R  M  A  T  < IX , "WATER*  ) 

0  02  3  3  1  FGRMATUX,  “FIELD") 

■'  o  2  9  3  2  FORMAT  (  IX, "FOREST "  > 

0  0  3  o  8  3  =■  0  R  H  h  T  <  1  ,  "  C  I  T  Y  "  ) 

0  0  3  1  8  4  FORMAT;;  IX,  “THIS  CARTOGRAPHIC  FEATURE  IS  HuT  SPEC  I  FI  ED  "  '• 

0032  338  I  Ft  LUOT  E 0  6 >  W  R I T  £ ( LUOT ,  880) 

v  0  3  3  8  8  0  F  0  R  fl  A  T  <  “  1  ■  > 

0 0  7  4  RETURN 

< 1  o  3  3  E  N  0 


